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Abstract 

Starting from the Schrodinger-equation of a composite system, we 
derive unified dynamics of a classical harmonic system coupled to an 
arbitrary quantized system. The classical subsystem is described by 
random phase-space coordinates entangled with the quantized vari- 
ables of the complementary subsystem. Our semiclassical equation 
is true in a sense that its predictions are identical to those of the 
fully quantized composite dynamics. This exact method applies to 
a broad class of theories, including e.g. the relativistic quantum- 
electrodynamics and the electron- fonon dynamics. 
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Introduction. A consistent model for dynamical coupling between clas- 
sical and quantum systems would be attractive for various theoretical and 
practical reasons. We shall deliberately call such model semiclassical. It 
would, for instance, explain the "collapse" of wave function during quantum 
measurement as a result of dynamical interaction between the classical mea- 
suring apparatus and the measured quantum system Furthermore, par- 
ticle physics would be understood on the classical background of space-time 
0. A remarkable theoretical proposal of a canonical semiclassical dynamics 
H was motivated by more practical aims, related to physics of molecules. 
Unfortunately, in Ref. [|IJ, some variables of the classical subsystem acquire 
quantum uncertainties. The mean-field model is, on the contrary, un- 
able to explain the obvious statistical uncertainties arising in the classical 
subsystem's dynamics due to the backreaction of the quantum subsystem. 
The canonical formalism j3| may not preserve the positivity of state distri- 
bution function, cf. Ref. M. As yet, all speculations have failed to obtain a 
consistent semiclassical model. 

For the time being, we are not able to model the backreaction of a quan- 
tum system onto a generic classical one. Here, nevertheless, harmonic oscil- 
lator is a favorable exception: its quantum and classical dynamics map onto 
each other exactly f5|. Classical phase-space coordinates appear as labels of 
coherent quantum states pi which form an overcomplet basis for the oscilla- 
tor's Hilbert-space. One inclines therefore to guess that classical description 
could be used for a system of quantum harmonic oscillators (harmonic sys- 
tem, in short) even if it interacts dynamically with another quantum system, 
provided the interaction is linear or at most quadratic in the harmonic vari- 
ables. In our Letter we show that, using coherent state representation, this 
is indeed the case. 

We derive semiclassical equation for the wave function of a composite 
system containing a harmonic subsystem. This subsystem is described in 
classical stochastic terms while it interacts with the quantized dynamics of the 
complementary subsystem. Our semiclassical equation provides, by quantum 
and stochastic averages, the true quantum expectation values for the fully 
quantized composite system operators. 

The model. We consider a quantum system consisting of two interact- 
ing subsystems. One of them is assumed to be harmonic. Its interaction 
with the other subsystem should be linear in the harmonic variables. This 
class includes many systems of primary physics interest like, e.g., quantum- 
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electrodynamics (hence quantum-optics as well), the electron-fonon interac- 
tion, the spin-boson system 0, or the so-called Caldeira-Leggett model [§]. 
We shall use the terminology of quantum-electrodynamics or quantum-optics 
whenever concrete terms help explanation. Accordingly, we consider the fol- 
lowing form of Hamiltonian: 

n n 

where the Planck's constant h has been set to 1. The absorption 
and creation operators for the harmonic system (which will be called field). 
The a> n 's are the field mode frequencies. Hq stands for the Hamiltonian of 
the charged particles (which we call atomic system, for simplicity's). The 
atomic operators J n coupled to the field modes will be called currents. In 
interaction picture, we get the following Schrodinger-equation for the state 
vector \l/ of our composite quantum system: 

= "«E [UtWn +ji{t)a n ) % (2) 

n 

where j n {t) = exp(iuj n t)J(t) is the "rotated" version of the interaction picture 
current J n (t) for each field mode. 

Instead of the abstract state vector we introduce a wave function in 
coherent state representation for the harmonic subsystem, i.e. for the field: 

tp(a,a*) = (a,a*\ty, (3) 

where \a, a*) is equal to the direct product of the coherent states |a n ,a*) 
of all field oscillators (n — 1,2,.. .), respectively. The scalar product on the 
RHS is to be taken on the factor Hilbert-space of the field. (The author has 
failed to find a more suitable compact notation.) Accordingly, ijj(a, a*) is 
wave function of the field coordinates a, a* while, on the other hand, it is 
abstract state vector in the atomic Hilbert-space. 

Conditional quantum state. We are going to give a special interpretation 
to the wave function (3). Let us consider the norm function which is known 
as the Q-fwactioia in quantum optics 0: 

w(a,a*) = ||^(a,a*)|| 2 . (4) 

Formally, it is a probability distribution over the classical phase-space of the 
field. From now on, let us choose this way to interpret it! Consequently, 
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the field will be considered classical and stochastic, described by random 
phase-space coordinates a, a*. 

Given an arbitrary function F(a,a*), its stochastic mean will reproduce 
the exact expectation value of the corresponding antinormally ordered oper- 
ator :F(a, a)): ant i as calculated in the original quantum state, i.e.: 

tf:F(a,a% nti V = [ F(a,a>(a, «*) II —• (5) 

This relation follows from standard literature on coherent states ||. It is 
straightforward to inspect that a more general relation holds true in cases 
where F(a, a*) is not necessarily an (a, a*)-dependent number but an arbi- 
trary atomic operator. Then we can write: 

¥:F(a,a% n u*= U\a,a*)F(a,a*)^a,a*)\{^. (6) 

Invoking the Eq. (4) and its probabilistic interpretation suggested above, 
it is reasonable to introduce the notion of conditional quantum expectation 
value of a given atomic operator F: 

= ^(a,a*)Fjj(a,a*) 
{ |^(«,a*)|| 2 ' 1 ' 

which, in fact, is the quantum expectation value of F at fixed state (a, a*) of 
the classical field. Lastly but most importantly, we shall interpret the wave 
function ip(a, a*) as conditional quantum state of the atomic subsystem at 
fixed classical state of the field. That this way of interpretation is consistent 
we see by applying Eqs. (4) and (7) on the RHS of Eq. (6): 

& :F(a, a% nti * = / (F(a, a*)) a ^w{^ a*) — • (8) 

In summary, the atom's state is described by state vector while the field's 
state is characterized by random phase-space coordinates (a, a*). The state 
of the composite system is given by the unnormalized conditional atomic 
state ip(a,a*) whose norm yields the phase-space distribution of the field. 
The original quantum expectation values of entangled atomic and normal 
ordered field operators can exactly be obtained by stochastic averaging of 
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the corresponding conditional atomic quantum expectation values over all 
classical fields. 

Equation of motion. In coherent state representation the following oper- 
ator correspondences should be used (cf. Ref. 0): 

4* <- a *if>(a, a*), a n V <- ( ^ + ij,(a, a*). (9) 

Hence the Schrodinger-equation (2) leads to the following equation of motion 
for the wave function (3): 

Ma, a*) = -i (a*j(t) + |jt(t)) ^(«, a*) - ij\t) d ^f\ (10) 

where we have suppressed the notation of indices n as well as of summation 
over them. 

From the standard theory of coherent states it follows that the conditional 
state vector (3) can always be written in the following form: 

if)(a,a*) = e" |Q|2/ V(a*). (11) 

Since the so-called Bargmann-state </?(«*) [0 does not depend on a the 
equation of motion (10) becomes simpler in the Bargmann-representation: 

&(a*) = -ia*j(t)tp t (a*) - yt(f)M^. (12) 

Discussion. To interpret the interaction between the the quantized atomic 
system and the classical field, respectively, we rewrite the equation of motion 
(10) [or (12), equivalently] for the conditional density operator defined by 
p(a, a*) = ip(a, a*)^(a, a*). Using the Eq. (11), too, one obtains: 
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a*j(t) + aj\t),p t ] -ij\t)^ + i^j{t). (13) 



The first term on the RHS is a Hamiltonian contribution to the quantum 
state dynamics as if the classical field were an external parameter. The 
derivative terms differ by operator orderings from the terms proposed in 
Ref. |§. The difference leads to consistency. These derivative terms repre- 
sent a specific kind of backreaction of quantum fluctuations, controlling the 
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stochastic spread of the classical system's variables. To see that, let us trace 
both sides of the equation over the atomic Hilbert-space. We obtain the 
following drift equation for the phase-space density (4): 

d 

w t (a, a*) = -i-Q^ (0' t (*))a,«* tt, t( a J a *)) + cx - • ( 14 ) 

Hence, the phase-space coordinates of each member in the random ensemble 
obey to deterministic classical Hamiltonian equations of motion 

a = -<0'(*)>«*«. > «* = *O' f 00U* . ( 15 ) 

where quantized atomic currents appear via their conditional quantum mean 
values. These equations coincide formally with the mean-field equations 
(cf.,e.g., in Refs. [0, |J). In mean- field theory the classical phase-space vari- 
ables take sharp values. It does not include the backreaction of the fluc- 
tuations in the quantum subsystem. In our theory, however, the influence 
of those fluctuations is reflected by the stochasticity of classical phase-space 
coordinates. 

Summary, outlook. We derived dynamic coupling between a classical 
harmonic system and a generic quantum system. Our semiclassical equation 
yields, via quantum and stochastic averages, all quantum expectation values 
of the full quantized dynamics of the composite system. The scope of this 
semiclassical theory is not restricted to the non-relativistic quantum mechan- 
ics. It applies to a broad class of theories, including first of all the relativistic 
quantum-electrodynamics, the electron-fonon interaction in solid states, the 
spin-boson problem etc. 

Admittedly, our claim that exact quantum results may be computed from 
stochastic averages would provoque concerns. Especially, the existence of a 
stochastic model for the radiation field interacting with quantized currents 
raises conceptual issues such as, e.g., a possible conflict with Bell's theorem 
which forbids the existence of local stochastic models. Here we do not in- 
tend to discusse the corresponding loopholes. Some measurement theoretical 



implications have already been discussed elsewhere |TT[ . 

We recall the fundamental issues mentioned in the introduction as main 
motivations. The collapse of wave function can certainly be described if we 
apply our equations to a suitably chosen model of measurement situation. 
Most interestingly, a natural collapse dynamics appears whenever we let our 
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quantum system to interact with the classical electromagnetic and/or grav- 
itational fields, described of course in the mathematical framework of the 
present Letter. Finally, we emphasize the practical value of the possibility 
that a great deal of interacting quantum systems becomes exactly treatable 
in semiclassical framework which is in many respect simpler than the original 
quantum one. 

Appendix: Principles of numeric simulation. We outline the principles of 
an efficient algorithm to simulate the solutions of the semiclassical equations 
of motion. Let us assume that we know the (Bargmann-) conditional state 
<fit( a *) a ^ some initial time t. Then the initial phase-space distribution of the 
field is 

w t (a 1 a*) = e-^ 2 \\i Pt ( a *)\\ 2 . (16) 

For the first numeric step, we generate a long random sequence of complex 
phase-space coordinates a(l), a(2), . . . , a(N) according to the distribution 
(16). Meanwhile, we keep the steps \a n (k + 1) — a n (k)\ small for all n and 
k. (Especially in case of many field modes, the Monte-Carlo "importance 
sampling" [0 seems useful to generate such random sequences.) We store 
the sequence {a(k); k = 1, N}. Then we calculate tp(k) = ip t (a*(k)) for each 
k and store the sequence {(p(k);k = 1,N}, too. These two sequences will 
represent the state of our system numerically. Their time evolution will be 
deterministic. Let us calculate their values for time t + e where the time 
increment e is small. From the equation of motion (15) we read out the 
updated phase-space coordinates: 

a n (k) -> a n (k) - (* = 1, 2, . . . , N; n = 1, 2, . . .). (17) 



The numeric approximation of the equation of motion (12) leads to the fol- 
lowing update of the conditional state vectors: 

„(*) **)-«£<«(*)*.(«>**) " fcgdM jg^jlgl) ' (18) 

for n = 1, 2, . . .. The upper sign is to be taken for k = 1, 2, . . . , iV — 1; we 
take the lower one for k = N. Repeated applications of the basic update 
cycle (17,18) yield the state of the system for any later time. The quantum 
expectation value of any (antinormally ordered) composite operator can be 
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calculated numerically. To calculate quantum expectation values, we approx- 
imate the RHS of the analytic relation (7) numerically as follows: 

-n^aUnu^--^ (19) 

i.e., it is approximately equal to the ensemble average of the corresponding 
conditional quantum expectation values. The relation becomes identity in 
the limit iV — > oo. 

The numeric derivation 

dyt(a*(k)) w y> t (q*(fc±l))-y> t (a*(fc)) 
(9a* a„(A; ± 1) — a n (fc) 

is, as usual, an awkward step of the algorithm. We have initially chosen small 
enough increments \a n (k + 1) — a n (k)\ to assure a good quality of numeric 
derivation. However, the "chain" {a(k);k = 1,2, ...,N} may, after many 
updates, becomes distorted is such a way that its increments are not uniform 
and small anymore. Then, a special formatting procedure must be applied 
to both sequences {a(k)} and {ip(k)} to restore uniform small increments 
before the updating steps continue. 

The author is indebted to J. Butterfield for his illuminating remark on 
conditional quantum states. This work was supported by the grants OTKA 
No. 1822/1991 and T016047. 
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